The Case for an Accelerating Universe from Supernovae 

To Appear as an Invited Review for PASP 



Adam G. RiessQ 
ABSTRACT 

The unexpected faintness of high-redshift Type la supernovae (SNe la), as measured by two 
teams, has been interpreted as evidence that the expansion of the Universe is accelerating. We 
review the current challenges to this interpretation and seek to answer whether the cosmological 
implications are compelling. We discuss future observations of SNe la which could offer 
extraordinary evidence to test acceleration. 
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1. Introduction 

Two teams have presented observational evidence from high-redshift type la supernovae (SNe la) that 
the expansion of the Universe is accelerating, propelled by vacuum energy (Riess et al. 1998; Perlmutter 
et al. 1999). The primary evidence for this hypothesis is the faintncss of distant SNe la relative to their 
expected brightness in a decelerating Universe. The question we propose to answer in this review is whether 
the observations of distant supernovae compel us to conclude that the expansion is accelerating. 



2. Supernova Measurements 
2.1. Past Work 

Supernovae (SNe) have a long history of employment in the quest to measure Hubble's constant (see 
Branch 1998 for review) and currently provide a column of support for a strong consensus that H o =60 to 
70 km s" 1 Mpc . The history of utilizing supernovae to measure the time evolution of the expansion rate 
is far briefer. All initial proposals for using high-redshift SNe to constrain global deceleration recognized 
the necessity of an optical space telescope to collect the data. Wagoner (1977) envisioned application 
of Baade's method or the the Expanding Photosphere method (Kirshner & Kwan 1974) to measure the 
angular diameter distance of Type I (hydrogen-deficient) and Type II (hydrogen-rich) supernovae at 
z = 0.3. Colgate (1979) demonstrated even greater prescience, suggesting that SNe I at z = 1 could be 
used as standard candles for "determining the cosmological constant with greater accuracy than other 
standard candles." Further thoughts by Tammann (1979) included the necessity to account for the redshift 
of the observed light using spectrophotometry of SNe in the ultraviolet (i.e., if -corrections) , host galaxy 
extinction, and time dilation of the light curves. 

Yet even before the launch of the Hubble Space Telescope (HST), a persistent two-year ground-based 
effort by a Danish group was rewarded by the discovery of their first (and only reported) high-redshift SN 
la, SN 1988U at z = 0.31 (Norgarrd-Nielsen et al. 1989) as well as modest bounds on the deceleration 
parameter, —0.6 < qo < 2.5. This team employed modern image processing techniques to scale the 
brightness and resolution of images of high-redshift clusters to match previous images and looked for 
supernovae in the difference frames. Unfortunately, the project's low discovery rate coupled with the 
dispersion of SNe la when treated as perfect standard candles (~0.5 mag) suggested that the determination 
of the deceleration parameter would require a scientific lifetime. 

However, great progress was made by the Supernova Cosmology Project (SCP) in the detection rate 
of high-redshift SNe la by employing large-format CCDs, large-aperture telescopes, and more sophisticated 
image-analysis techniques (Perlmutter et al. 1995). These advances led to the detection of seven SNe la at 
z w 0.4 between 1992 and 1994, yielding a confidence region that suggested of a flat, A = universe but 
with a large range of uncertainty (Perlmutter et al. 1997). 

The High-z Supernova Search Team (HZT) joined the hunt for high-redshift SNe la with their discovery 
of SN 1995K at z = 0.48 (Schmidt et al. 1998). Both teams made rapid improvements in their ability to 
discover ever greater numbers of SNe la at still larger redshifts. One-time redshift record holders included 
SN 1997ap at z = 0.83 (SCP; Perlmutter et al. 1998), SN 1997ck at z = 0.97 (HZT; Garnavich et al. 1998a), 
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SN 1998eq at z = 1.20 (SCP; Aldering et al. 1998) and SN 1999fv at z=1.23 (HZT; Tonry ct al. 1999). 
Before either teams' claimed samples of high-redshift SNe la were large enough to detect the acceleration 
signal, both teams found the data to be inconsistent with a Universe closed by matter (Garnavich et al. 
1998a; Perlmutter et al. 1998). 



Record redshift Type la Supernovae 



SN la 


z 


discoverers 


discovered 


SN 1988U 


0.31 


N0rgaard-Nielsen-Niclsen et al; IAUC 4641 


8/1988 


SN 1992bi 


0.46 


SCP; IAUC 5652 


4/1992 


SN 1995K 


0.48 


HZT; IAUC 6160 


3/1995 


SN 1995at 


0.66 


SCP; IAUC 6270 


12/1995 


SN 1996cl 


0.83 


SCP; IAUC 6621 


3/1996 


SN 1997ap 


0.83 


SCP; IAUC 6596 


3/1997 


SN 1997ck* 


0.97 


HZT; IAUC 6646 


4/1997 


SN 1998eq 


1.20 


SCP; IAUC 7046 


10/1998 


SN 1999fv 


1.23 


HZT; IAUC 7312 


11/1999 



* Poor spectrum of SN but light curve consistent with SN la at measured z of host 

These observational feats were preceded by increased understanding and ability to make use of SNe la 
observations to constrain the cosmological parameters. Empirical correlations between SN la light curve 
shapes and peak luminosity improved the precision of distance estimates beyond the standard candle model 
(Phillips 1993; Hamuy et al. 1995; Riess, Press, Kirshncr 1995; Perlmutter et al. 1995; Tripp 1997; Saha et 
al. 1999; Parodi et al. 2000). Studies of SN la colors provided the means to distinguish supernovae which 
were reddened by dust from those which were intrinsically red (Riess, Press, Kirshner 1996; Riess et al. 
1998; Phillips et al. 1998; Tripp & Branch 1999). Goobar & Perlmutter (1995) showed that measurements 
of SNe la at different redshift intervals could break degeneracies between and Oa. Additional work on 
cross-filter ^-corrections provided the ability to accurately transform the observations of high-redshift SNe 
la to the rest-frame (Kim, Goobar, & Perlmutter 1996). 



2.2. Observations 

In order to resolve whether the results from high-redshift SNe la are compelling it is important to 
review how the measurements were obtained. 

Both the SCP and the HZT detected their samples of high-redshift SNe la by using large-format CCDs 
at telescopes with large apertures (most commonly the CTIO 4-m Blanco telescope). Using two sets of 
deep images in the R or /-band spaced across a lunation, the "template" images are subtracted from the 
second-epoch images and automated software searches for sources in the difference images whose intensity 
surpasses a specific threshold. The observations taken in pairs are spaced over a small time interval of a 
few minutes to eliminate moving transients. Human inspectors "filter" the automated results in an effort 
to maximize the likelihood that candidates are supernovae (Schmidt et al. 1998; Perlmutter et al. 1999). 
Because the resources available for collecting spectral identifications for all the candidates are insufficient, 
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candidates which appear most likely to be SNe la are given priority. The factors which are favored in the 
human selection criteria include candidates which are separated from the host galaxy nuclei and those 
that show good contrast with the host galaxy (especially those with little or no apparent host). The 
signal-to-noise ratios of the identifying spectra vary greatly (Perlmutter et al. 1995, 1998; Riess et al. 1998) 
but have improved with the availability of the Keck Telescope and an increased emphasis on the search for 
clues of SN la evolution. Most of the spectral identifications were made by visual comparison to template 
spectra of nearby SNe la. More recently automated cross-correlation techniques have been employed (Riess 
et al. 1997). Approximately half of the SN la redshifts were determined from narrow emission lines or Ca 
H and K absorption; the rest were derived from the broad supernova features. Although formal statistics 
are not currently available, we are aware that the SCP candidates have yielded a greater fraction of SNe la 
than the HZT candidates, a significant number of which turn out to be SNe II. Taken together, about half 
of the two teams' candidates are revealed to be SNe la, with the rest classified as SNe II, AGNs, flare stars, 
or unclassified objects. SNe la in the desired redshift range are monitored photometrically in two colors by 
observatories around the world using red-sensitive passbands; -ffSThas monitored the light curves of about 
two dozen of these objects. 

The initial search template images are eventually replaced by deeper template images obtained with 
good seeing and taken one or more years after discovery when the SNe la have faded by five to eight 
magnitudes. More recently, repeated searching of the same fields has provided deep template images 
before the SN explosion. After the SN images are bias corrected and flat-fielded they are geometrically 
aligned, and the resolution and intensity are scaled to match those of the templates. After subtracting 
the host galaxy light the SN magnitudes are measured. The SCP uses aperture photometry; the HZT 
fits point-spread-functions. Uncertainties are determined synthetically by the injection of artificial SNe of 
known brightness. The SCP uses standard passbands and Landolt (1982) standards of comparison while 
the HZT uses a custom passband system which is transformed to the Landolt scale (Schmidt et al. 1998). 
Custom cross-band if-corrections are calculated using spectrophotometry of nearby SNe la whose colors 
are reddened to match the high-redshift objects. 

The distances are measured by fitting empirical families of light curves to the flux observations of 
individual supernovae. The measured distances are derived from the luminosity distance, 



where C and T are the SN's intrinsic luminosity and observed flux, respectively and Dl is in Mpc. 
Alternately, a logarithmic luminosity distance (i.e., the distance modulus) is used: 



where M is the SN's absolute magnitude and m is the observed magnitude in a given passband. 

Three different light-curve fitting methods have been used to measure the distances, each of which 
determines the shape of the best-fitting light curve to identify the individual luminosity of the SN la. The 
HZT has used both the multi-color light-curve shape method (MLCS; Riess, Press, & Kirshner 1996; Riess 
et al. 1998) and a template fitting approach based on the parameter Amis(i?) (Phillips 1993; Phillips et 
al. 1999), while the SCP uses the "stretch method" (Perlmutter et al. 1995; 1999). 

Nearby SNe la provide both the measure of the Hubble flow and the means to calibrate the relationship 
between light-curve shape and luminosity. The SCP uses ~20 nearby SNe la in the Hubble flow from the 




(1) 



./i = m - M = 51ogD L + 25 



(2) 



- 5 - 



Calan/Tololo Survey (Hamuy et al. 1996) while the HZT adds to this set an equal number of SNe from the 
CfA Sample (Riess et al. 1999). 

Dust extinction is handled somewhat differently by the two teams. The HZT measures the extinction 
from the B — V reddening (i.e., the color excess) and then combines this measurement in a Bayesian 
formalism using a prior host galaxy extinction distribution calculated by Hatano, Branch, & Deaton (1998). 
This treatment assumes that extinction makes supernovae appear dimmer (farther), never brighter (closer). 
The SCP uses a number of different approaches (including the HZT approach), but favors making no 
individual extinction corrections and discarding outliers. Figure 1 shows a single Hubble diagram made 
with the data from both teams (Riess et al. 1998; Perlmutter et al. 1999). 

The measured distances are then compared to those expected for their redshifts as a function of the 
cosmological parameters Qm, fiA, and H n : 

D L = cH^(l + z) \fl k \- 1/2 smn{\n k \ 1/2 f dz[(l + zf{\ + Q M z) - z{2 + z)fi A r 1/2 }, (3) 



where fife = 1 — fijvf — fiA, and sinn is sinh for fife > and sin for fife < (Carroll, Press, & Turner 1992). 
The likelihoods for cosmological parameters are determined by minimizing the \ 2 statistic between the 
measured and predicted distances (Riess et al. 1998; Perlmutter et al. 1999). Determination of fiju and 
fiA are independent of the value of H n or the absolute magnitude calibration of SNe la. The confidence 
intervals in the fiM-fiA plane determined by both teams are shown in Figure 2. In this plane, acceleration 
is defined by the region where the deceleration parameter, qo, is negative, 

<Zo = ^-fiA<0. (4) 



2.3. Results 



The visual impression from Figure 2 is that the supernova observations favor the parameter space 
containing a positive cosmological constant and an accelerating Universe with high statistical confidence. 
The confidence regions determined by the two teams are in remarkably good agreement. Both teams claim 
that these results do not result from chance with more than 99% confidence (Riess et al. 1998; Perlmutter 
et al. 1999). However, the specific confidence of these cosmological conclusions depends on which parameter 
space intervals one considers to be equally likely, a priori. This point is addressed by Drell, Loredo, & 
Wasserman (2000) who suggest that models with fiA=0 and fiA ^ could be considered equally probable, 
a priori. While Riess et al. (1998) and Perlmutter et al. (1999) considered the probability that fiA > 
with a flat prior in a linear fiA space, an alternative is to use a flat prior in a logarithmic fiA space. 

Another useful way to quantify the SN la constraints has been given by Perlmutter et al. (1999) as 
0.8fiM — 0.6fiA = —0.2 ± 0.1, a result which applies equally well to the Riess et al. (1998) data. 

A more illuminating way to quantify the evidence for an accelerating Universe is to consider how the 
SN la distances depart from decelerating or "coasting" models. The average high-redshift SN la is 0.19 mag 
dimmer or ^10% farther than expected for a Universe with no cosmological constant and negligible matter 
(fiA = 0, fijw = 0). Of course it's apparent that the Universe has more than negligible matter and the 
current consensus from the mass, light, X-ray emission, numbers, and motions of clusters of galaxies is that 
fi M « 0.3 (Carlberg et al. 1996; Bahcall, Fan, & Cen 1997; Lin et al. 1996; Strauss & Willick 1995). The 
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high-redshift SNe la from both teams are 0.28 mag dimmer or 14% farther than expected in a Universe 
with this much matter and no cosmological constant. The statistical uncertainty of these values is 0.08 and 
0.06 mag (or 4% and 3% in distance) for the HZT and SCP, respectively. The observed dispersion of the 
high-redshift SNe la around the best-fit cosmology is 0.21 mag for the HZT and 0.36 mag for the SCP. A 
frequentist would consider the accelerating Universe to be statistically likely at the 3cr to 4ct level (i.e., > 
99%), while the Bayseian likelihood would depend on a statement of the natural space and scale for Q^. 

Simply put, high-redshift SNe la are ~0.25 mag fainter than expected in our Universe with its presumed 
mass density but without a cosmological constant (or which is not accelerating). The statistical confidence 
that SNe la are fainter than expected is high enough to accept that it does not result from chance and 
additional SNe la continue to support this conclusion (Schmidt 2000). Rather, this result is only challenged 
by systematic uncertainties not reflected in the variance of high-redshift SN la distance measurements. In 
the following sections we review the challenges to the cosmological interpretation of the SN la observations 
and consider whether the evidence compels us to believe that the Universe is accelerating. 



3. Challenges and Tests of the Accelerating Universe 
3.1. Evolution 

Could SNe la at z = 0.5, a look-back time of ~ 5 Gyr, be intrinsically fainter than nearby SNe la by 
25%? For the purpose of using SNe la as distance indicators near and far, we are concerned only with 
an evolution which changes the luminosity of a SN la for a fixed light-curve shape. Evolution is a major 
obstacle to the measurement of cosmological parameters, having plagued workers who tried to infer the 
global deceleration rate from brightest cluster galaxies in the 1970s (Sandage & Hardy 1973). We will 
consider both the theoretical and empirical indications for SN la evolution. 

Theoretical understanding of SNe la provides reasons to believe that evolution is not a challenge to the 
accelerating Universe. SNe la are events which occur on stellar scales, not galaxian scales, and therefore 
should be less subject to the known evolution of stellar populations. However, our inability to conclusively 
identify the progenitor systems (see Livio 2000 for a review) and our lack of a complete theoretical model 
(see Leibundgut 2000 for a review) means we cannot rely exclusively on theory to rule out the critical 
degree of evolution. 

Nevertheless, theoretical calculations can provide some insight into this question. Hoflich, Wheeler, & 
Thielemann (1998) have calculated models of spectra of SNe la with solar and one-third solar metallicities 
and have found little difference between the spectral energy distribution over the wavelengths where the 
SNe la have been observed (see Figure 3). In principle, changes in the age and hence initial mass of the 
progenitor star at high redshift could yield white dwarfs of varying carbon-to-oxygen (C/O) ratio. It is 
currently difficult to assess if such a variation could produce significant evolution as these calculations lack 
the necessary precision (Dominguez ct al. 1999; von Hippel, Bothun, & Schommer 1997). Umeda et al. 
(1999) suggest that a change in the C/O ratio is the source of the inhomogencity in SN la luminosity, but 
they conclude that calibration of the luminosity via light-curve shape relations effectively inoculates the 
cosmological measurements to an evolution in the C/O ratio. 

To date, answers to the question of whether SNe la evolve have been sought from empirical evidence. 
In the nearby sample, SNe la are observed in a wide range of host-galaxy morphologies including ellipticals, 
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post-starburst galaxies (e.g., SN 1972E in NGC 5253), low surface brightness galaxies (e.g., SN 1995ak in 
IC 1844), irregulars (e.g., SN 1937C in IC 4182), SOs (e.g., SN 1995D in NGC 2962), and early to late-type 
spirals (van den Bergh 1994; Cappellaro et al. 1997; Hamuy et al. 1996; Riess et al. 1999). The range of 
metallicity, stellar age, and interstellar environments probed by the nearby hosts is much greater than the 
mean evolution in these properties for individual galaxies between z = and z — 0.5. Some variation in 
the observed characteristics of SNe la with host morphology has been seen in the nearby sample (Hamuy 
et al. 1996; Branch, Romanishin, & Baron 1996) Yet after correction for the light-curve shape/luminosity 
relationship and extinction, the observed residuals from the Hubble flow do not correlate with host galaxy 
morphology (see Figure 4). In addition, Hubble flow residuals show no correlation with the projected 
distance from the host center (see §3.6). This evidence suggests that SN la distance estimates are insensitive 
to variations in the supernova progenitor environment and is the strongest argument against significant 
evolution to z = 0.5 (Schmidt et al. 1998). However, this evidence is still circumstantial as we cannot 
be sure that the local environments of the SN la progenitors arc similar to the average environments of 
the hosts. Future studies which probe the local regions of nearby SNe la should be able to explore their 
variance. 

The other empirical test of evolution has been to compare the observed characteristics of low-redshift 
and high-rcdshift SNe la. The assumption of this test is that a luminosity evolution of ~25% would be 
accompanied by other visibly altered characteristics of the explosion. Here too our lack of firm theoretical 
footing makes it difficult to gauge the correspondence between any evolution in distance-independent 
quantities and luminosity. Therefore we must conservatively demand that observations of all observables of 
distant SNe la be statistically consistent with the nearby sample. 



3.1.1. Spectra 

Comparisons of high-quality spectra between nearby and high-redshift SN la, such as those seen in 
Figure 5, have revealed remarkable similarity (Riess et al. 1998; Perlmutter et al. 1998, 1999; Filippenko 
et al. 2000). The spectral energy distribution is sensitive to the atmospheric conditions of the supernova 
(i.e., temperature, abundances, and ejecta velocities). Even primitive modeling indicates that it would 
be difficult to retain the primary features of the SN la spectrum while altering the luminosity by about 
20% to 30%. Further, comparisons of temporal sequences of spectra reveal no apparent differences as the 
photosphere recedes in mass (Filippenko et al. 2000), indicating that the superficial similarities persist at 
deeper layers. However, among the variety of nearby SNe la are objects which are both ~25% fainter than 
the average and also display very typical spectral features (e.g., SN 1992A, see Figure 5). Therefore, the 
existing spectral test alone is not sufficient to check for this degree of evolution. 

While spectral similarity between nearby and distant SNe la provides no indication of evolution, a 
lack of any spectral peculiarities among high-redshift SNe la could signal some changes at high redshift. 
Li et al. (2000) find from the most unbiased survey of nearby SNe la to date that ~ 20% of SNe la are 
spectroscopically similar to the overluminous SN 1991T. SN 1991T-like objects show weak Ca II, Si II, and 
S II, but prominent features of Fc III (Filippenko 1997), and close cousins, such as SN 1999aa, have similar 
characteristics with the exception that Ca II absorption is more normal. Monte Carlo simulations of the 
search criteria used by the SCP and the HZT team performed by Li, Filippenko, & Riess (2000) indicate 
that such overluminous objects should comprise approximately 25% of high-redshift SNe la (with some 
uncertainty due to a possible link between such objects and circumstellar dust). To date, neither team has 
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reported the existence of a single SN 1991T-like object among ~ 100 high-rcdshift objects. 

It is certain that the low signal-to-noise ratio of the spectra of high-redshift SNe la, coupled with 
the rcdshifting of spectral features out of the observable window makes it more difficult to identify 
individuals from this peculiar class. In addition, the spectroscopic peculiarities of SN 1991T-like objects 
are only apparent close to maximum light (or earlier), and some high-redshift SNe la may not have been 
observed early enough to identify their spectral peculiarities. This same effect may also explain why the 
Calan/Tololo survey of 29 SNe la yielded no SN 1991T-like objects (Hamuy et al. 1996). If, however, these 
observational biases are not to blame, the absence of SN 1991T-like SNe la at high redshift could result 
from an evolution of the population of progenitor systems (see Livio 2000 for a review) or a subtle difference 
in selection criteria (see §3.6). If true, this type of evolution may yield important clues which help identify 
the progenitor systems (Ruiz-Lapuente & Canal 1998), but it is unlikely to affect the measurement of the 
cosmological parameters since spectroscopically normal SNe la at low and high redshift have been used to 
derive the cosmological constraints. 



3.1.2. Broad-band 

The distributions of light-curve shapes for nearby and distant SNe la are statistically consistent 
(Riess et al. 1998; Perlmutter et al. 1999). Such consistency appears to extend to infrared light curves of 
high-redshift SNe la which show the characteristic second maximum of typical, low-redshift SNe la (Riess 
et al. 2000). 

An analysis by Drell et al. (2000) indicates that different light-curve fitting methods may not be 
statistically consistent and that the apparent differences may be a function of the light curve shape. 
However, these conclusions are highly sensitive to estimates of the correlated distance uncertainties between 
different fitting methods and these correlated uncertainties are difficult to estimate. 

Evolutionary changes in the model temperature and hence the thermal output of the explosion could 
be detected from the colors of pre-nebular supernovae. The most significant analysis of B — V colors, 
performed by Perlmutter et al. (1999), demonstrated consistency between low and high-redshift SNe la 
at maximum light. Likewise, Riess et al. (2000) found that the B — I colors of a SN la at z = 0.5 were 
consistent with those of nearby SNe la. However, Falco et al. 1999 (see also McLeod et al. 1999) suggested 
that the B — V colors of high-rcdshift SNe la from the HZT may be excessively blue, a conclusion which 
cannot be rejected by the B — I color measurements by Riess et al. (2000). More data are needed to confirm 
or refute this possibility. If true, this could indicate either evolution or the existence of a halo of Milky 
Way dust which would redden the observed wavelengths of nearby SNe la more than redshifted objects. 
This latter possibility has been suggested by recent Milky Way dust maps of Schlcgcl, Finkbcincr, & Davis 
(1999) in contrast to the previous maps of Burstein & Hciles (1982), but it would augment rather than fully 
explain the faintness of distant SNe la. 

The risetime (i.e., the time interval between explosion and maximum light) is sensitive to the ejecta 
opacity and the distribution of 56 Ni . The risetime of nearby SNe la (Riess et al. 1999b) and the 
high-rcdshift SNe la of the SCP (Goldhaber 1998; Groom 1998) were initially strongly discrepant (Riess et 
al. 1999c). However, a reanalysis of the SCP high-rcdshift data by Aldering, Knop, & Nugent (2000) finds 
the high-redshift risetime to be longer and much more uncertain than indicated by Groom (1998). The 
remaining difference could be no more than a ~2.0er chance occurrence. More early photometry of distant 
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SNe la is needed to increase the significance of this test of evolution. 

Evolution is arguably the most serious challenge to the cosmological interpretation of high-redshift 
SNe la. Further studies, currently underway, seek to compare the host galaxy morphologies and luminosity 
versus light-curve shape relations for nearby and distant SNe la. The results reviewed in this section do 
not appear to provide any clear evidence of evolution. However, absence of evidence is not necessarily 
evidence of absence. The paucity of high signal-to-noise ratio observations of high-redshift SNe la and the 
current lack of a comprehensive theoretical model or a well-understood progenitor system keeps the embers 
of skepticism aglow. 



3.2. Dust 

Consideration of a non-cosmological explanation for the dimming of distant supernova light must 
invariably turn to a famous pitfall of optical astronomy: extinction. Trouble has often followed dust in 
astronomy, a point first appreciated by Trumpler (1930) when analyzing the spatial distribution of Galactic 
stars. 



3.2.1. Ordinary Dust 

An additional ^25% opacity of visual light by dust in the light paths of distant supernovae would be 
sufficient to nullify the measurement of the accelerating Universe. Both teams currently measure SN la 
colors to correct for the ordinary kind of interstellar extinction which reddens light. Galactic extinction 
maps from Burstein & Hcilcs (1982) and Schlegel, Finkbeiner, & Davis (1998) were used by the HZT 
and the SCP, respectively, to correct individual SNe la for Milky Way extinction. Such corrections were 
typically less than 0.1 mag due to the high galactic latitudes of the SNe la. Even a previously unknown 
halo of Galactic dust would dim the restframe light of nearby SNe la more than highly redshifted SNe la 
and would therefore not explain the cosmological indications. 

Measurements of B — V colors have been used by both teams to test for and remove host galaxy 
extinction (Riess et al. 1998; Perlmutter et al. 1999; see Figure 6). Totani & Kobayashi (1999) have 
suggested that the remaining uncertainty in the mean measured B — V color excess (a — 0.02 mag; 
Perlmutter et al. 1999), when multipled by reddening ratios of 3 to 4 to determine the optical opacity, 
may be too large to discriminate between open and lambda-dominated cosmologies with high confidence. 
However, such concern seems unwarranted as this uncertainty remains 3 to 4 times smaller than the size of 
the cosmological effect of an accelerating Universe. 

Another measurable effect of the critical amount of mean interstellar extinction is that it would 
introduce more dispersion in the distance measurements than is currently observed. A random line of sight 
into a host galaxy will intersect a nonuniform amount of extinction. Hatano, Branch, & Deaton (1998) 
have calculated the expected distribution of extinction along random lines of sight into host galaxies. A 
mean, uncorrected extinction of 0.25 mag would induce twice the distance dispersion observed by the HZT 
(Riess et al. 1998). In addition, high-redshift surveys are biased towards finding SNe la which have even 
less extinction than would be expected from the distributions of Hatano et al. (1998). 
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A more powerful way to search for reddening by dust is to observe high-rcdshift SNe la over a large 
wavelength span: from the optical to the infrared. Infrared color excesses would be more than twice as 
large as E B -v for ordinary dust. A set of such observations for SN 1999Q (z = 0.46) disfavor A v = 0.25 
mag of dust with Galactic-type reddening at high confidence (Ricss et al. 2000), but more SNe la need to 
be observed in the near-infrared to strengthen this conclusion. 



3.2.2. Grey Dust 

More pernicious than ordinary dust is "grey" dust which could leave little or no imprint on the spectral 
energy distribution of SNe la. Perfectly grey dust is only a theoretical construct, but dust which is greyer 
than Galactic- type dust (i.e., larger reddening ratios) does exist (Mathis 1990) and could challenge the 
cosmological interpretation of high-redshift SNe la. 

Grey dust can be made with large spherical dust grains or elongated "whiskers." Past studies 
of whiskers (Aguirre 1999a; Rana 1979, 1980) indicate that they would distort the cosmic microwave 
background (CMB), an effect which has not been seen. Like non-grey extinction, grey interstellar extinction 
does not provide an acceptable explanation for the dimness of SNe because the inherent variations in the 
opacity along random lines of sight would induce more distance dispersion than is observed (Riess et al. 
1998). 

Grey intergalactic extinction could affect measurements of the deceleration parameter (Eigenson 1949) 
without tell-tale dispersion or reddening. Indeed, observations of neither SNe la nor other astrophysical 
objects rule out a 30% opacity by large semi-spherical dust grains (Aguirre 1999b). Aguirre (1999a, b) 
has shown that a uniformly distributed component of intergalactic grey dust with a mass density of 
Qdust ~ 5 x 10~ 5 and graphite grains of size > 0.1 ^im could explain the faintness of high-z SNe la without 
detectable reddening and without overproducing the currently unresolved portion of the far-infrared (far-IR) 
background. However, this physical model of dust would provide some reddening which can readily be 
detected with observations in the optical and infrared. Measurements by Ricss et al. (2000; sex; Figure 7) of 
Eb-i for a single high-redshift SN la disfavor a 30% visual opcaity of grey dust at the ~ 2.5cr confidence 
level, but more observations are needed to strengthen this conclusion. Additional studies of the faint far-IR 
sources seen with SCUBA may soon provide definitive constraints on the unresolved component of the 
far-IR background and the viability of extragalactic grey dust. 



3.3. Gravitational Lensing 

The inhomogeneous distribution of matter in the Universe typically deamplifies and very rarely 
amplifies the observed brightness of distant SNe la compared to the average. (Note that the mean observed 
brightness must equal the unamplified value expected in a perfectly smooth Universe.) 

The size of the typical deamplification is a function of the SN la redshift, the mass density of the 
Universe and the fraction of dark matter locked into compact objects. This effect has been quantified by 
a wide range of techniques (Kantowski, Vaughan, & Branch 1995; Frieman 1997; Wambsganss et al. 1997; 
Holz & Wald 1998; Kantowski 1998; Metcalf 1999; Barber 2000). The effect of weak lensing on the observed 
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distribution of luminosities of SNe la at z — 1 and z — 0.5 can be seen in Figure 8. In the most relevant 
regime for the current SNe la at z = 0.5 (i.e., CIm ~ 0.3, and mostly diffuse dark matter) the typical 
dcamplification is ^2%, much smaller than the cosmological effect. An extreme case (i.e., Ojy « 0.5, all 
matter in point masses) could deamplify the median SN la at z = 0.5 by 5% (Holz 1998), but this model 
is unlikely to be correct and the effect is still not large enough to negate the cosmological interpretation of 
high-redshift SNe la. Perlmutter et al. (1999) considered lensing by up to Dm = 0.25 in compact material 
in the determination of their confidence intervals. They found little impact on the likelihood of a positive 
cosmological constant (see Figure 9). Wang (2000) suggests that by flux-averaging (i.e., binning the SNe la 
distances by redshift) one can reduce the bias due to weak lensing. In the future, any bias due to weak 
lensing will naturally vanish as the sample sizes become larger and the mean observed luminosity more 
robust. It is interesting and potentially useful to note that the observed distribution of SN la distances (see 
Figure 8) can in principle be used to determine the fraction of gravitating matter contained in compact 
objects (Seljak & Holz 1999; Metcalf & Silk 1999). 



3.4. Measurement Biases 

In this section we consider if biases in the measurement process of high-redshift SNe la could mimic the 
evidence for an accelerating Universe. An exhaustive list of such biases has been considered by Hogg (2000) 
and Hogg & Turner (1998). Here we discuss how these biases may apply to the supernova measurements. 

The observational challenge is to measure the distance to high-redshift SNe la which are 6 to 7 
magnitudes fainter and have lower signal-to-noise ratio than those which delineate the Hubble flow. 
Differences in the way low and high-redshift SNe la are observed must not introduce biases in their 
distance measures at more than the few percent level. Indeed, Hogg (2000) has noted that the proximity of 
high-redshift SNe la to any reasonable world model is a testament to the feasibility of measuring distances 
across such a large range. However, because the goal of these observations is precision cosmology (and not 
simply to demonstrate the dynamic range of useful photometry) our scrutiny must be greater. 

Charge transfer inefficiency (CTI) and detector non-linearities can cause faint objects to appear fainter. 
However, ground-based observations of high-redshift SN la are limited by the bright sky, a regime in which 
these effects are widely found to be negligible. For space-based observatories such as the Hubble Space 
Telescope, CTI is far more troublesome but quite correctable (Whitmore, Heyer, & Casertano 1999). In 
addition, only a subset of the high-redshift SNe la have been measured with HST and the cosmological 
conclusions do not depend on the inclusion of these objects. 

High precision, flux-conserving algorithms have been developed to properly subtract images of the 
host galaxy from images with SN light (Alard & Lupton 1998). Correctly employed, these methods reduce 
any biases in the measurement of the SN brightness to less than a few percent. Tests for measurement 
biases and estimates of uncertainty are performed by both teams by the injection of artificial SNe into the 
observed images. 

Hogg & Turner (1998) discuss a bias towards higher observed fluxes which naturally occurs when 
measuring the brightness at discovery of low signal-to-noise ratio sources. This bias results from the 
preferential selection of faint sources on the bright side of the Poisson distribution of photon statistics. 
Follow-up observations of the source would not incur this bias. This effect would have little impact on the 
supernova distances measured by the HZT and SCP because the light curves are dominated by observations 



- 12 - 



made after discovery. In addition, the direction of this effect is opposite to the signal of an accelerating 
Universe. However, this effect may become more important for SNe la found at z > 1 for which the 
discovery observation may provide one of the most significant measurements. 



3.5. Selection Biases 

Do the HZT and SCP preferentially select faint SNe la at high redshift? Because we have already 
considered evolution in §3.1, here we are only concerned with the characteristics of a high-redshift sample 
which is drawn from the same population as the nearby sample. In so doing, we must also consider if the 
nearby sample is a fair representation of that population. 

As an example, consider the set of nearby SNe la which appear fainter than expected for their redshift 
in the bottom panel of Figure 1. Presumably these objects appear dim due to the intrinsic random scatter 
of SNe la. If, however, these SNe la had a characteristic in common which, in addition, favored their 
discovery at high-redshift, a bias would result. To date, no such characteristic has been identified and the 
observed dispersion of nearby SNe la is consistent with their measurement errors. 

Howell, Wang, & Wheeler (2000) found a difference between the projected distances from the hosts' 
centers for the nearby and distant SNe la (see Figure 10). Many of the nearby SNe la were found in the 
photographic Calan/Tololo survey in which saturated galaxy cores masked SNe near their hosts' centers 
(Shaw 1979; Hamuy & Pinto 1999). The result is that distant SNe la are more centrally located than the 
nearby sample. However, in an analysis of 44 nearby SNe la, Rless et al. (1999) found no dependence of 
the distance measurement on the projected distance from the host center, so this selection effect appears to 
have no bearing on the cosmological use of SNe la (see Figure 11). 

Malmquist bias (Malmquist 1924; 1936) can shift the mean distance too close in a magnitude-limited 
survey of SNe la. This effect seems to contrast with the cosmological dimming perceived in an accelerating 
Universe. However, if the nearby sample were more afflicted by Malmquist bias than the distant sample, 
this bias could mimic an accelerating Universe. Because the intrinsic scatter of SN la distances is low 
(< 0.15 mag), Malmquist bias, which scales with the square of the dispersion (see Mihalas & Binney 1981 
for a derivation), is small for SNe la. Perlmutter et al. (1999) made analytic calculations of Malmquist bias 
arising from the intrinsic dispersion of SNe la (assumed to be 0.17 mag) and the SCP search incompleteness 
(determined empirically) to estimate that the net bias between the samples is no more than 0.03 mag. 
[Perlmutter et al. (1999) notes that the net bias may actually be closer to zero due to a compensating bias 
against the selection of light curves which are "fast" for their luminosity and therefore spend less time 
above the detection limit.] Ricss et al. (1998) used a Monte Carlo exercise to simulate the selection of SNe 
la near and far. Inputs to this exercise included the time interval between successive search epochs, limiting 
magnitudes, observed light-curve shapes, and the distribution of SN la luminosities. They report a net 
bias of less than 0.01 mag. These results indicate that the net Malmquist bias has negligible impact on the 
cosmological conclusions. 



3.6. Alternative Cosmological Models 
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The conclusions drawn from high-rcdshift SNe la are predicated on a model with two free parameters, 
Qm and f^A, and a Friedmann- Robertson- Walker cosmology. In the absence of a sound fundamental 
motivation for w Qmi alternate and more general descriptions of an energy density with negative pressure 
have been suggested (Caldwell, Dave, & Steinhardt 1998). These phenomenological or "quintessence" 
models invoke a decaying scalar field rolling down a potential as the source of today's acceleration (Wang 
ct al. 2000). A distinction of these models from a cosmological constant is that w, the ratio of pressure 
to energy density, is between -1 and 0, whereas w is exactly -1 for a cosmological constant. For feasible 
quintessence models, w, the equation-of-state parameter, varies slowly with time and can be approximated 
today by a constant equal to 



where a is the scale factor and Qq is the energy density of the vacuum component. The current acceleration 
for these models (assuming only two significant energy components today, Qm and Qq) is 



and is generally less than for a cosmological constant (all other parameters fixed). Inspection of equation 
(2) reveals that the Universe is accelerating if q is negative (d(t ) is positive), requiring that w < — |, 
independent of the value of Qm ■ 

Can we determine if the expansion is accelerating in a quintessence model? The SN la data from 
Pcrlmutter et al. (1999) and Riess et al. (1998) already provide meaningful constraints on w (Garnavich 
ct al. 1998b; Perlmutter et al. 1999). Increasing the value of w from -1 (i.e., for a cosmological constant) 
reduces the acceleration provided by a fixed value of Qq, but larger values of Qq are needed to retain an 
acceptable fit to the data. Graphically, increasing w from -1 rotates the error ellipses in Figure 2 to favor 
lower values of Qm and greater values of Qq. As seen from equation (6), the line separating an accelerating 
and decelerating Universe rotates in the same direction (always anchored at Qm = and Qq = 0), providing 
no gain on an acceptable region of parameter space which is not accelerating. The nearest intersection 
between a non-accelerating region of parameter space and one which is preferred by the data remains when 
Q m « 1 and Qq << 1. However, values near Qm = and Qq = are poor fits to the data independent 
of the value of w. 

Perhaps the simplest way to understand why the SN la data favor an accelerating Universe is to 
consider an FRW cosmology with Q M = and no vacuum energy. This empty Universe must be neither 
accelerating nor decelerating but simply coasting. The fact that the high-redshift SNe la are systematically 
farther for their redshift than expected in this cosmology means that the distance between low and 
high-redshift SNe la (where redshift is a surrogate for time) grew faster than expected for a Universe which 
has been coasting on today's Hubble expansion. This implies that the Universe has been accelerating. 

An alternate cosmological explanation to acceleration has been posited by Goodwin et al. (1999) and 
Tomita (2000). They suggest that the supernova data are also consistent with a decrease in the Hubble 
expansion by 10-20% beyond z = 0.1 (300/i _1 Mpc). The distance at which the Hubble expansion dips 
would correspond to the approximate radius of the "local" underdensity in which we live. Although a few 
peculiar flow surveys support bulk motions on scales up to half this size (Lauer & Postman 1994; Hudson 
et al. 1999), most recent surveys do not (Dale et al. 1999; Courteau et al. 2000; Colless et al. 1999; Riess 
et al. 1999; see Willick 2000 for a review). However, the biggest problem with such a commodious, local 
underdensity is its great improbability. Power spectra demonstrate (Watkins & Feldman 1995; Feldman 
& Watkins 1998) that the density of the Universe is extremely homogeneous on this scale, and finding 




(5) 



q = -a(t a )a(t )/a 2 (t a ) = -(Q M + Qq(1 + 3io)) 



(6) 
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ourselves in the midst of such a vacuous location would be virtually anti-Copcrnican. Using cold dark 
matter power spectra constrained by CMB observations and large scale structure, Shi and Turner (1998) 
and Wang, Spergel, and Turner (1998) expect 0.5% to 1.5% variations in the Hubble constant on 300/i _1 
Mpc scales, a factor of 20 times smaller than required in the local void model. By filling in the Hubble 
diagram of SNe la at 0.1 < z < 0.2 it would be possible to test this model. 

Outside the FRW cosmologies the SN la data can have significantly different interpretations. For 
example, in steady-state cosmologies, SN redshifts do not come from expansion, but rather through 
"tired-light" processes. However, the SN la data exhibit the time dilation effect expected in an expanding 
Universe, implying that the tired-light hypothesis is incorrect (Leibundgut et al. 1996; Goldhaber et al. 
1996; Riess et al. 1997; but see Narlikar & Arp 1997). In the quasi-steady state cosmology, the SN la 
data lead to modifications of the model, such as matter creation during periodic expansion phases (Hoyle, 
Burbidgc, & Narlikar 2000). A detailed consideration of how to interpret the SN la data in non-FRW 
cosmologies is beyond the scope of this review, but is thoroughly addressed by Hoyle, Burbidge, & Narlikar 
(2000). Alternative theories of gravity such as modified Newtonian dynamics models (MOND; Milgrom 
1983, 1998; McGaugh & de Blok 1998a, b) could also modify the interpretation of the observations of 
high-redshift SNe la. 



4. Conclusion 

After reviewing the cosmological interpretation of SN la observations and the current challenges to the 
analysis of the data, we can now offer an answer to the question initially posed: do the observations of 
distant supernovae compel us to conclude that the expansion of the Universe is accelerating? 

With full consideration of the evidence, we conclude that an accelerating Universe remains the most 
likely interpretation of the data because the alternatives, individually, appear less likely. However, the 
quantity and quality of the SN la evidence alone is not yet sufficient to compel belief in an accelerating 
Universe. The primary sources of reasonable doubt are evolution and extinction, as discussed above. 
Although the types of studies also described above could potentially yield evidence that either of these 
non-cosmological contaminants is significant, the current absence of such evidence does not suffice as 
definitive evidence of their absence. Our current inability to identify the progenitors of SNe la and to 
formulate a self-consistent model of their explosions exacerbates such doubts. Even optimists would 
acknowledge that neither of these theoretical challenges is likely to be met in the near future. 

Fortunately there are at least two routes to obtain compelling evidence to accept (or refute) the 
accelerating Universe, one of which employs the use of SNe la at even greater redshifts. 



5. Epilogue 

5.1. The Era of Deceleration 

If the Universe is accelerating, it is a rather recent phenomenon likely commencing between z « 0.4 
and 1. Before this time the Universe was more compact and the pull of matter dominated the push of 
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vacuum energy in the equation of motion. As a result, the Universe at z > 1 must be decelerating. This 
cosmological signature should be readily apparent by extending the Hubble diagram of SNe la to z « 1.2. 
By this redshift SNe la in an accelerating Universe will cease to diverge in distance from an equally massive 
Universe without vacuum energy. Alternatively, if a monotonically increasing, systematic effect is the 
source of the excessive faintness of high-redshift SNe la, the measured distances of SNe la at z > 1 will 
continue to diverge from a cosmology without vacuum energy and in addition would diverge from the 
cosmological model inferred from SNe la at z = 0.5 (see Figure 12). Complex parameterizations of evolution 
or extinction selected to match both the accelerating and decelerating epochs of expansion would require a 
near conspiracy of fine tuning and are highly doubtful. 

Efforts are already underway to find and measure SNe la at z > 1. Gilliland, Nugent, & Phillips 
(1999) used a subsequent epoch of the Hubble Deep Field to detect two SNe, one (SN 1997ff) with a 
photometric redshift of z = 1.32. The elliptical host of SN 1997ff suggests that this object is of Type la, 
but the observations are insufficient to provide a useful distance estimate. The SCP reported the discovery 
of SN 1998ef at z = 1.20 (Aldering et al. 1998) and follow-up observations with the HST will provide a 
useful distance estimate (Aldering 2000). The HZT recently reported the discovery of four SNe la at z > 1 
including SN 1999fv at z = 1.23 (Tonry et al. 1999). From this growing sample will likely come the means 
to search for the epoch of deceleration. 



5.2. Cosmic Complements 

We previously sought to determine if the observations of SNe la alone require an accelerating Universe. 
Now we will briefly consider the cosmological constraints provided by other astrophysical phenomena. A 
thorough discussion of these constraints is beyond the scope of this review but can be found elsewhere 
(Tyson & Turner 1999; Roos & Harun-or-Rashid 2000) 

Current measurements of the cosmic microwave background (CMB) power spectrum indicate that 
the sum total of energy densities is within 10% of unity. This result is seen from the BOOMERANG 
(Melchiorri et al. 1999; Lange et al. 2000; dc Bcrnardis et al. 2000), the TACO (Miller et al. 1999) and 
MAXIMA experiments (Hanany et al. 2000; Balbi et al. 2000) and from a compilation of all other CMB 
measurements (Tegmark & Zaldarriaga 2000). In addition, estimates of Qm from the mass, light, X-ray 
emission, numbers, and motions of clusters of galaxies converge around 0.2 to 0.3 (Carlberg et al. 1996; 
Bahcall, Fan, & Cen 1997; Lin et al. 1996; Strauss & Willick 1995). These two pieces of information alone 
indicate a significant contribution by vacuum energy, sufficient to produce an accelerating Universe (see 
Figure 13). Additional constraints from observations of the Lyman-alpha forest, cluster evolution, double 
radio galaxies, and statistics of gravitational lenses have been used to tighten these conclusions (Roos & 
Harun-or-Rashid 2000; Turner 1999; Eiscnstcin, Hu, & Tegmark 1999; Lineweaver 1998). 

Although no single cosmological observation yields a conclusive census of the energy densities in 
the Universe, the combined constraints from multiple experiments is providing strong bounds on the 
cosmological parameters. Each individual experiment has unique sources of systematic uncertainty. By 
combining the results of many experiments it should be possible to negate their impact on the determination 
of the cosmological parameters. 



I wish to thank Mario Livio, Alex Filippenko, and Robert Kirshncr for helpful discussions and 
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Figure Captions: 

Figure 1: Hubble diagrams of SNe la from Perlmutter et al. (1999; SCP) and Riess et al. (1998; HZT). 
Overplotted are three cosmologies: "low" and "high" VLm with J7a = and the best fit for a flat cosmology, 
57m = 0.3, f^A = 0.7. The bottom panel shows the difference between data and models from the Hm = 0.3, 
= prediction. The average difference between the data and the Q,m = 0.3, J7a = prediction is 0.28 
mag. 

Figure 2: SNe la joint confidence intervals for (Um,^a) from Perlmutter et al. (1999; SCP) and Riess 
et al. (1998; HZT). Regions representing specific cosmological scenarios are indicated. 

Figure 3: Type la supernova model spectral energy distributions for solar and ^ metallicities. 
Superimposed are the transmission functions for standard passbands; from left to right is U, B, V, and R. 
From Hoflich, Wheeler & Thielemann (1998) 

Figure 4: The nearby Hubble diagram of SNe la in different host galaxy morphologies. In the top 
panel the SNe la are treated as standard candles and in the bottom panel distances are determined with 
the MLCS method (Riess, Press, & Kirshner 1996; Riess et al. 1998). After MLCS corrections are made, 
the distance Hubble flow residuals are independent of host galaxy morphology. 

Figure 5: Spectra of 4 nearby and 1 high-redshift SN la at the same phase (Riess et al. 1998). Within 
the observed variance of SN la spectral features, the spectra of high-rcdshift SNe la are indistinguishable 
from the low-redshift SNe la. The spectra of the low-redshift SNe la were resampled and convolved with 
Gaussian noise to match the quality of the spectrum of SN 1998ai. 

Figure 6: Comparison of color excess, Eb-v, f° r nearby SNe la from the Calan/Tololo survey (Hamuy 
et al. 1996a, b) and high-redshift SNe la from Perlmutter et al. (SCP; 1999). The color excesses at high and 
low redshift are consistent indicating negligible extinction and/or no evidence for color evolution. 

Figure 7: The color evolution, B — I, and the color excess, Eb-i, of a h igh-redshift SN la, SN 1999Q 
(Z = 0.46), compared to the custom MLCS template curve with no dust and enough dust (of either 
Galactic- type or greyer) to nullify the cosmological constant. The smaller error bars are from photometry 
noise; the larger error bars include all sources of uncertainty such as intrinsic dispersion of SN la B — I 
color, X-corrections, and photometry zeropoints. The data for SN 1999Q are consistent with no reddening 
by dust, moderately inconsistent with ^4y=0.3 mag of grey dust (i.e., graphite dust with minimum size 
> 0.1 //m; Aguirre 1999a, b) and Ay=0.3 mag of Galactic-type dust. From Riess et al. (2000). 

Figure 8: Probability distribution P(p) for supernova apparent brightness \i normalized to /i = 1 for a 
filled beam (i.e., a homogeneous Universe). The vertical lines are at the empty-beam value. "Galaxies" are 
treated as isothermal spheres and truncated at a radius of 380 kpc; "compact objects" are point masses. 
From Holz (1998.) 

Figure 9: Cosmological constraints from Perlmutter et al. (SCP; 1999) for three weak lensing scenarios. 
Fit C assumes a filled beam, case K assumes an empty beam, and fit L is a model with weak lensing by up 
to VLm = 0.25 in compact objects. 

Figure 10: The projected distances from the host centers of nearby SNe la discovered photographically 
and with CCDs and high-redshift SNe la discovered with CCDs (Howell, Wang, & Wheeler 2000). 

Figure 11: MLCS distance residuals from the Hubble flow for nearby SNe la versus their projected 
distance from their host centers (Riess et al. 1999a). 



Figure 12: The Hubble diagram of SNe la (see Figure 1) and the effect of a systematic error which 
grows linearly with redshift (e.g., evolution or grey extinction). 

Figure 13: Cosmological constraints from SNe la, CMB, and matter (Turner 1999). 
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